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Using numerical modeling we study emergence of structure and structure-related nonlinear con- 
duction properties in the self-assembled nanoparticle films. Particularly, we show how different 
nanoparticle networks emerge within assembly processes with molecular bio-recognition binding. 
We then simulate the charge transport under voltage bias via single-electron tunnelings through the 
junctions between nanoparticles on such type of networks. We show how the regular nanoparticle ar- 
ray and topologically inhomogeneous nanonetworks affect the charge transport. We find long-range 
correlations in the time series of charge fluctuation at individual nanoparticles and of flow along 
the junctions within the network. These correlations explain the occurrence of a large nonlinear- 
ity in the simulated and experimentally measured current-voltage characteristics and non-Gaussian 
fluctuations of the current at the electrode. 

PACS numbers: 

I. INTRODUCTION 

Functional soft materials composed of many nano-particles often appear in complex structure both as 2-dimensional 
films til, i2j and 3-dimensional super-crystal structures t^, iSj. These structures emerge in various processes of self- 
assembly, most of which utilize liquid degrees of freedom and nonequilibrium methods [1, 0] ■ The gold nanoparticles 
are often used due to the property that their surfaces can be functionalized in an appropriate manner. The functional- 
ization of the nanoparticle surfaces has effects on their interactions [8] . In the case of nanoaprticle films on substrate 
assembled by the liquid evaporation Q the thiol-passivated nanoparticles where used, where the attached chains 
determine the minimum distance between two neighbouring nanoparticles. Another way to induce desired interaction 
between two nanoparticles was designed recently with attaching the pieces of biologically compatible strands of DNA 
molecules to different nanoparticles The DNA- mediated interaction then utilizes the biological recognition and 
binding energy of the compatible DNA base pairs. Depending on the number of other factors, which we will also dis- 
cuss later, this programmable self-assembly may lead to random as well as regular super-crystal structures, as shown 
in recent experiments in Refs. 0, Somewhat similar (but in a different potential range), is the bio- recognition 
assembly which utilizes biologically compatible protein molecules, for instance the antigen-antibody recognition p^ . 
which can be used in biosensors. 

The nanoparticle assemblies exhibit a wide range of physical properties which are not found in the bulk materials 
and at the level of a single nanoparticle § . Moreover, the emergent physical properties of the nanoparticle assemblies 
can often be related to their spatial structure HJBiSllH- For instance, new optical and magnetic properties have been 
found in super-lattices made of nanocrystal B y sp eciffic functionalization of their surfaces gold nanoparticles 

develop directionally dependent exchange interaction [12| . Assuming that such nanoparticles can be arranged on nodes 
of a sparse scale- free graph, numerical simulations predict new features of the magnetic response of such a structure 
[l^ . [l4 | , compared to the classical memory materials [isj . Similarly, the films of gold nanoparticles under voltage 
bias conduct current via single- electron tunneling processes (l6| due to their small capacitances and the Coulomb 
blocade effects In contrast to classical conduction, the single-electron conduction processes have numerous 

advantages, e.g., the absense of dissipation due to heathing, and the nonlinear current-voltage characteristics, which 
are important in the nanoelectronics. Experimentally a large nonlinearity of the current-voltage curve has been 
found in self-assembled nanoparticle films on substrates, which appears to be fairly correlated with the film structure 

The random structure of nanoparticle films on substrates can be adequately modelled by planar graphs (nano- 
networks) [H, |2l| . For the charge transport via single-electron processes, which we study here, the topological elements 
of these nano-networks are suitably selected as follows: the nodes represent nanoparticles, whereas edges (links) are 
the tunneling junctions between two nanoparticles. For the junction (link) to occur between two nanoparticles, the 
distance between them must be smaller than the tunneling radius of the electrons for that kind of nanoparticles [13] ■ 
The charge transport through the film under applied voltage can be described within the model of the capacitively- 
coupled nanoparticles on substrate [23|. In our previous work ^23] we generalized the model of capacitively-coupled 
nanoparticle assemblies for an arbitrary topology of the nano-networks 18, 23] and implemented a numerical model 
of single-electron tunnelings driven by the external voltage through such networks. Theoretical background and the 
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details of the numerical implementation is reviewed in Ref . psj . The simulation results revealed the topology-induced 
nonlinearity in the I{V) curve, in a good agreement with the experiments [T8| . 

In the present work we study the collective dynamic behavior which makes the physical basis for the observed 
I{V) nonlinearity, in particular, the long-range correlations in charge fluctuations monitored at each nanoparticle and 
inhomogeneous charge flow along the conducting paths through the sample and at the electrode. We determine the 
probability distributions of various quantities related to the conduction process and demonstrate how the topological 
and charge disorder affect the quantitative statistical properties and the I{V) dependences. 

The organization of the paper is as follows: In Section 2 we present the model of the particle self-assembly with 
bio-recognition bonding and show some of the typical emergent structures. In section 3 we summarize the theoretical 
background of the conduction with single-electron processes and its numerical implementation for an arbitrary network 
topology. The visualization of the conduction paths through two network types, the regular and inhomogeneous 
nanoparticle aggregates, as well as the statistical properties of the topological and charge flow are shown in Section 
4. In Section 5 we analyse the simulated time series of charge and current fluctuations and the related I(V) curves. 
Section 6 gives a short summary of the results. 



II. NANOPARTICLE NETWORKS EMERGING IN SELF-ASSEMBLY WITH BIO-RECOGNITION 



Different assembly processes make use of the liquid phase, that provides the nanoparticle diffusion and affects the 
interaction. The assembly with fast evaporation of the liquid P, Q induces additional constraints to the particle 
motion and interaction. The interaction itself between pairs of the nanoparticles moving in a liquid is taking part 
at short distances and depends on the type of functionalization of their surfaces, as mentioned above. In the case 
of bio-recognition binding, which we will discuss shortly here, the attached DNA induce the differences between two 
types of nanoparticles, i.e., type A and type B nanoparticles, each of which carry the biologically compatible strand of 
the same DNA molecule. Thus, the interaction between A-B type is attractive within some bonding distance d. The 
bonding distance depends on the lengths of the attached molecules and on the number of the compatible base-pairs 
on them ^^]. Therefore, it exceeds the sum of the particle radii dh > Ra + Rb- The effective interaction between A-B 



particles can be then modeled by the Lennard- Jones potential [24 1 



where r is the actual distance between particles and the potential minimum corresponds to the optimal bonding 
distance. In the present case, the parameters a and e depend, roughly speaking, on the lengths of the attached DNA 
chains and the strength of binding (number of the compatible base pairs) [rj . On the contrary, the particles which 
carry biologically incompatible DNA strands have no attractive interaction of this type. The repulsive part alone 
of the potential in Eq. ([T]) can describe the interaction of the equal- type particles A-A and B~B, occurring at short 
distances. 

The diffusion of such nanoparticles can be simulated via the Langevin equation with the random force originating 
from the integrated liquid degrees of freedom and the thermal fluctuations [25| . Apart form the molecular difference 
between attached DNA strands, generally, the particle radii Ra ^ Rb can be different and consequently two different 
diffusion coefficient enter the equation 

l'^r^ = -Vr. ^ V^-j{\r, - T^]) + F,^T, (2) 



with (Fi^rit)) = 0, {Fi^TFj^xit')) = GkBTi/iSijSt.t' , and Di = kBT/vi, i referring to either A or B particle type. A 
complete description of the model and the discussion of the role of different parameters in the self-assembly processes 
will be given elsewhere [l^. Here we restrict our discussion to the situation which is interesting for the assembly 
of the conducting nanoparticle films. Therefore, we consider particles of equal size Ra = Rb, and vary the relative 
concentration ca '■ cb of particles A and B carrying with different strands of the molecule. In this way various 
structures can be assembled, both regular and irregular, which in turn appear to have different physical properties 
(as discussed below). 

In Fig. [1] we show two representative structures which emerge in the process when temperature is kept low compared 
to the parameter e, and the total coverage of the surface by particles is fixed, but the relative concentrations of A and 
B particles varied. Generally, for low enough temperatures and equal concentrations the system has tendency towards 
a regular ordered structure with few defects (Fig. [1] right). Whereas, when the concentrations are in disballance, 
another type of structure with local order and large voids may appear, as shown in Fig.[Tl left. It should be stressed 
that, for a fixed set of parameters, many such structures may appear depending on the initial conditions. Each of 
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FIG. 1; Two examples of the emergent nanoparticle films obtained via bio-recognition binding between A (red) and B (yellow) 
particles for equal coverage and different relative concentrations ca/cb as follows 25% : 75% (left) and 50% : 50% (right). 



them corresponds to a local free energy minimum, and appear to be quite stable for the low temperatures [25l |. In the 
following we will discuss in detail the conduction properties of a regular nanoparticle structure and of a highly irregular 
nanoparticle structure on the substrate. In practice, the regular (triangular) nanoparticle arrays jlS] and irregular 
nanoparticle films Q are assembled using entirely different approaches. Our results shown above demonstrate that 
such diverse structures may appear in the same self-assembly process with bio-recognition when a single parameter 
(the relative concentration of differently functionalized nanoparticles) is suitably varied. 



III. SINGLE-ELECTRON TUNNELINGS IN NANOPARTICLE FILMS 



The capacitively coupled nanoparticle assembly on substrate under applied voltage conducts current with a cascade 
processes of single-electron tunnelings, which is enabled by the Coulomb blockade at junctions between the adjacent 
nanoparticles [itL [2^. Following a single electron tunneling through the junction between two particles a — > 6, the 
charging energy of the nanoparticle increases inversely proportional to its capacitance. This increases the barrier for 
further tunnelings along that junction, and consequently, no tunnelings occur until the voltage is increased enough to 
overcome the barrier. 

The theoretical background of the tunnelings within the nanoparticle arrays has been developed in [l^, . Recently 
we have generalized the formulation for the case of nano-networks with the tunneling junctions representing as the 
network edges [l^. It should be stressed that the links in the nano-network suitable for the tunneling processes 
might differ form the ones of the bio-recognition binding, shown in Fig. [1] Specifically, for the tunneling processes 
important are positions of the nanoparticles on the substrate and their mutual distances compared to the tunneling 
radius rt- The chemical bonds, which are shown in Fig. [TJ are to be understood as fixing the nanoparticle positions. 
Then the tunneling junction network may have the links which coincide with the chemical bonds only when the 
bonding distance dh vt, otherwise, the tunneling junction network is different. Generally, the nano-network edge 
representing the tunneling junction occurs between any two nanoparticles which are separated at distances r < rt. 
In addition, the difference between A and B particle types, which was important for the assembly, plays no role in 
the conduction (assuming the particles of the same size and chemical composition). In this way, a nano-network of 
tunneling junctions can be constructed whenever the positions of the nanoparticles on the substrate are known [2^ . 
In our numerical models the tunneling radius is a free parameter. Varying the tunneling radius, e.g., by changing the 
chemical composition and size of the nanoparticles, slightly different nano-networks of junctions can be constructed. 
For the numerical implementation of the tunneling process the nano-network of tunneling junctions is fixed by the 
adjacency matrix with the elemnts Aij = 1, when the tunneling between the nanoparticles i and j can occur, and 
Aij = 0, when the junction is absent, e.g., due to large distance dij > rt between them. 
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A. Electrostatic energy of a nanoparticle network 

Following the classical work [22j the electrostatic energy for a capacitively coupled nanoparticle array consisting 
of N nanoparticles and of a general topology given by the adjacency matrix {^^j} can be described as follows [23j . 
The capacitances are related to the structure of the underlying nano-network by = CA^j. The charge Qi on zth 
nanoparticle can be expressed as 

Q« = E - + E - (3) 

j M 

where $i is the potential of the ith nanoparticle in the assembly. The index /i G {+,—,g} is attributed to the 
electrodes (including the gate) and thus $^ is the potential at the ^th electrode and similarly is the charge and 
Ci_f^ the capacitance between the electrode /ith electrode and ith the nanoparticle. The system of equations ([3]) for 
i = 1, 2, • • ■ can be written in the matrix form Q ~ — 0^$'' and formally solved for # as 

$ = M-iQ + A/"^C^$'' (4) 
where now M is the capacitance matrix of the whole system (nanoparticles plus electrodes) with the elements: 

k M 

Using these definitions, the electrostatic energy of the whole system can be written as 

£;= igtM-iQ + Q-V^^' + Q^*^, ^ M-^Cf,<^>''. (6) 

where fj, S {+, —,g} stands for the electrodes and gate, V'^^* is the external voltage and Ci^^, Q^, $^ are described 
above. 

B. Numerical implementation and visualization of charge transport 

The system is driven by increasing the external voltage at one, say left, electrode. Voltage on the gate is kept fixed. 
Due to the voltage bias, the electrons first tunnel from the left electrode to the neighboring nanoparticles and then 
forward along the available junctions between the nanoparticles. Due to the overall voltage profile the tunnelings are 
predominantly to the right, towards the zero-voltage electrode as +V —> a — > 5 ^ • • • ^ 0. 

Following the tunneling of a single electron from particle a b, causing the change in the particle charges as 
Q'i = Qi ~^ ^ib — Sia, the energy change AE(a — > b) occurs in the whole array, due the long-range interactions. The 
energy change is computed using the expression in Eq. ([5]) as follows: 

AE{a -^b) = Vt-Va + i(A4-' + AI,-,' Af^V - ^C) + - VT' (7) 

and similarly bE(a ±) = ±\/j + \M'^^ originating from the tunneling between the electrode and the nanoparticle. 
Here we have defined the following quantity 

T4^EQ»^^-\ (8) 

i 

which is also suitable for the numerical implementation. Namely, the updated value of the variable Yc at node c 
following the tunneling a ~> b can be computed recursively as follows: 

a^b: l/' = y, + Af-/-Af-i, a ^ ± : V: = Y±M-^\ (9) 

The simulation takes the following steps: For the fixed voltage V the whole network is updated for tunnelings. For 
each junction the attempted tunneling from i ^ j induces the global energy change AEi ^ j, which is computed 
using the above formulas. Then the tunneling rate for i ^ j event at given voltage V is computed according to (23j 

r. -(v)-— ^^--■'•/^ (10) 
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where e stands for the electron charge, T is the temperature and Ri^j is the quantum resistance of the junction. The 
quantum resistance (or conductance) can be computed for small nanoparticles. e.g., by density functional theory [27j . 
For instance, for the metallic nanoparticles one can obtain within the transfer Hamiltonian that Rt = 7. — oi„i(i^„ — =^ , 

where Tq and Dq are the transmission coefficient and the density of states at each nanoparticle [I4I • We assume that 
the quantum resistancei?^^^ = Rt for all junctions. Note that in the approach described above, Rt appears as a 
parameter which sets the characteristic time scale of the process via r 1/RtCg (see later). In our continuous time 
implementation, out of all attempted tunnelings only those which tunneling times ti^j taken from the probability 
distribution p{ti^j) = pQexp{—Ti^jt) a ppe ar to be in small time window St are processed. Then the network is 
updated again (for more details see Ref. |23l|).Here we use the limit C/Cg <C 1 and Vg — 0. 

When the external voltage is increased, more electrons tunnel into the nanoparticle array and consequently local 
ratio between the charge and voltage is changed, providing the conditions for more forward tunnelings. The moving 
front of charges on the slab of regular triangular arrays of nanoparticles is shown in Fig. IIIIBI for a set of different 
voltages. The color intensity indicates the number of electrones at the nanoparticle. When the voltage is large enough 
V ^ Vt, the moving front reaches the zero- voltage electrode and the current through the sample can be detected for 
the first time. For the voltages V > Vr the current through the sample is nonzero for all times. The average current 
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FIG. 2: Visualization of the charge profile before tlie voltage thresliold is readied, 
through the junction z — > j is defined as the net number of electrons per time unit: 



,{V)=eiT,^j{V) 



(11) 



In the numerical experiment we can measure the current through each junction within the film and at the electrode. 
In the laboratory experiments the current is measured only at the electrode. The current increases with increased 
voltage. However, the curren-voltage dependence appears to be non-linear for a range of voltages above the threshold 
V > Vt: 



I{V) ^ B 



V -Vt 
Vt 



C 



;C>1 ; 



(12) 



The exponent C, was found to vary from 2 to 4 with the structure of the nanoparticle arrays, in particular C, — 2.25 
was measured in the case of regular triangular array of gold nanoparticles (20| . whereas C — 3.9 and larger was found 
in highly inhomogeneous films obtained by the evaporation methods As mentioned in the introduction, one 
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of the purposes of the present work is to study the charge and current fluctuations inside the sample in order to 
reveal the correlations in the process which contribute the observed nonlinearity in the current-voltage curves. The 
other aspect is to demonstrate the role of random charge disorder in the transport process. The charge disorder 
represents fractional charges which are locally present in the real samples, for instance, due to differences in the 
contacts of some nanoparticles with the substrate, and affect the charge-voltage relation at these nanoparticles. In 
the simulations we consider the regular triangular array of nanoparticles with quenched charge disorder < < 1 at 
ith nanoparticle taken from the uniform random distribution. The calculations are systematically compared with the 
case of nanoparticle film without charge disorder but with large topological inhomogeneity of the structure (Fig. lIVp. 



IV. TOPOLOGICAL AND CHARGE FLOW: CONDUCTING PATHS 



For voltages over the threshold Vr the conducting paths through the sample are established. By recording the current 
through all junctions in the network we can visualize the charge flow along these conducting paths, as shown in Fig. lIVI 
Darker color corresponds to larger number of electrons transported along each link. The Fig. IIVI also demonstrates 
the role of topological inhomogeneity, i.e., in the case of the nano- network in the right panel, as opposed to the 
regular nanoparticle array. In the inhomogeneous nano-network, the areas without nanoparticles make a restriction 
to the charge transport, thus the number of conducting paths merge into a smaller number in the bottleneck regions. 
Draining the electrons along such paths thus affects a large regions behand. This also indicates that certain nodes 
on the paths play more important role in the charge flow than the others. The quantitative measure of the node 
importance is their topological betweenncss centrality B, which is defined as the number of shortest paths between 
all pairs of nodes on the graph that go through that node. In the case of the inhomogeneous film network in Fig. 
Iivr right) , the topological betweenness-centrality is given by a broad distribution with a power-law tail, shown in Fig. 
mjleft). In the regular nano-array, on the contrary, the topological betweenness-centrality of all nodes is similar (not 
shown), which induces a less intricate pattern of conducting paths. In Fig. lIVH eft) we show the conducting paths in 
the case of a triangular array in the presence of charge disorder. As already mentioned above, the presence of a random 
fractional charge at nanoparticle may destroy the local balance between the voltage and charge, leading to blocking of 
certain conducting paths through that nanoparticle and consequently reduced charge transport through that particle. 
Such areas with reduced charge flow appear as bright spots in Fig. lIVT left) . The distributions of the actual charge 
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FIG. 3: Visualization of the conduction paths through the regular triangular array with charge disorder at nodes (left) and 
through the topologically inhomogeneous network without charge disorder (right). For a fixed voltage V > Vt the cumulative 
charge flow along each fink is recorded and represented by the color intensity. 



flow (dynamical betweenness-centrality) for both networks is shown in Fig. [ll^right). Again the differences between 
two topologies play the role. In the inhomogeneous network NNETl, the distribution of the charge flow along the 
links on this network resembles the distribution of topological betweenness (one should stress that the shortest paths 
between the electrodes play the role in the conduction"). The simulation results can be approximately fitted with the 
g-exponential form [28|] P{F) - [1 - (1 - q){F/Fo)]^/'^^~'^'^ with q = 1.05. In the case of the triangular array with 
charge disorder, the distribution away from the peak area indicates that certain links in the network curry less charge 



7 




FIG. 4: Distribution of the topological flow B in the inhomogeneous nanonetwork NNETl (left) and distributions of charge 
flow F along the links in NNETl and in the triangular nanoparticle array with charge disorder (right). 



V. NON-GAUSSIAN CURRENT & CHARGE FLUCTUATIONS 



In our simulations the current recorded at the right electrode for the applied voltages above the threshold V > Vt 
can be compared with the experimental measurements in the nanoparticle films l2Clj|. The threshold Vr depends 
on the network size (distance between the electrodes) and its structure. 



A. Current— Voltage Curves 



In Fig.[5l[left) the current at the electrode is shown against reduced voltage {V — Vt) /Vt obtained in the simulations 
on different nano- net works. It is clear that for a wide range of voltages the dependences I{V) are nonlinear with 
C > 1 according to Eq. (fT2| for all 2-dimensional nano- networks. We have also shown the simulations for the 1- 
dimensional array of nanoparticles for the comparison. For the electron tunnelings through the 1-dimensional chain 
of nanoparticles the exponent remains C — 1 (top curve in Fig. [5]). However, in the 2-dimensional samples many 
conducting paths may open between the electrodes. Flow along these paths is a dynamical process with long-range 
correlations. The contributions of flow along different paths builds up the cumulative current which is measured at 
the electrode. Therefore, different geometries have different patterns of the conducting paths (as also shown in Fig. 
lIVp . leading to different nonlinear effects in the current-voltage dependences. Generally, the regular structures have 
smaller exponent, for instance, we find C ~ 2.25 for the triangular array with charge disorder, in very good agreement 
with the experiment [20|. A bit larger exponent is found in hexagonal nano-network without charge C, ~ 2.8, whereas 
in the inhomogeneous nanonetwork NNETl we have C, w 3.9, within the numerical error. The exponents in this range 
are found in other inhomogeneous nanoparticle films both in the simulations and in the experiment p^ . 

It is interesting that such cooperative effects of different paths occurs in a limited range of voltages V > Vr, 
approximately one order of magnitude on the reduced voltage scale (see Fig. [5]). In this range of voltages, the draining 
of the electrons along a conducting path is intermittent and induces an avalanche of flows in the connected area of 
the network. However, when the voltage is large enough, a permanent flow along each paths sets in, making them 
chain-like, and consequently the I{V) curve exhibits the crossover towards the linear dependence. 



B. Time-series Analysis 



In complex dynamical systems, time series of fluctuating variables contain valuable information about the correla- 
tions in the dynamical process [29|. In the case of the single-electron transport through the nanoparticle films, we 
consider temporal fluctuations of the current, which is measured at the electrode, and also the fluctuations of the 
charges (number of tunnelings) between all nanoparticles inside the sample. 
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FIG. 5: (left) Current plotted against reduced voltage for the topologically inhomogeneous nano-network NNETl (bottom 
curve) and for several regular nanoparticle arrays, as indicated: triangular array with charge disorder, and hexagonal and 
1-dimensional arrays, (right) Distributions of the increments of charge-fluctuations P{AI) plotted against A7 for the regular 
triangular array with charge disorder and for the topologically disordered NNETl. Fit are g-Gaussian according to Eq. (|13|l. 



For the current, which is measured at the electrode, we consider the time-series of the current {/(t)}, given in 
Fig. [6l and the increments of the current {AI{t)}, which is a measure of the acceleration in the electron transport 
along the conducting paths. Each point is averaged over time window twiN = 100 time steps (updates of the whole 
network). In Fig.[5l[right) we have shown the distributions of the increments P{AI) obtained for a fixed voltage above 
the threshold, for the triangular array with charge disorder and for the topologically disordered nanoparticle array 
NNETl. The distributions are fitted with the q-Gaussian [28l.[30t 



P(A/) = S(l-(l-g)(^ 



2\ 1/1-9 



(13) 



with different q values for different geometries of the nanoparticle arrays: q — 1.14 for the geometrically regular 
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FIG. 6; Current fluctuations time series for triangular nanoparticle array with charge disorder and for NNETl (left) and the 
corresponding power spectra (right). 



triangular array with charge disorder, while q — 1.85 for the topologically inhomogeneous NNETl. Similarly, the 
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analysis of the time series of the current itself, {I{t)}, exhibits different correlations in these two geometries. In Fig. [6] 
the power-spectra of the time series are shown. The spectrum of the topologically inhomogeneous NNETl shows the 
long-range correlations as S{i') ^ 1/v dX high frequencies, whereas the spectrum flattens at high frequencies in the 
case of the triangular charge-disordered array. At low frequencies, however, weak correlations are observed in both 
networks. 
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FIG. 7: Time series representing the number of tunnelings and their corresponding power spectra in (left) the triangular array 
with charge disorder and (right) in the topologically inhomogeneous NNETl (data from Ref. H^). 

The correlations in the current fluctuations at the electrode, i.e., at the last layer of the nanoparticles connected to 
the electrode, are caused by the tunneling processes deeper in the sample. The time-series of the number of tunnelings 
{Nxit)} at all nanoparticles in in the whole network and their power spectra are shown in Fig. [71 The power spectrum 
of the number of tunnelings exhibits long-range correlations for high frequencies, according to 

S{v) ^ v-^ ■ (14) 

The exponent (/> = 1 (within numerical error bars) for the topologically disordered nano- network NNETl, suggesting 
strong cooperation between the conducting paths. Whereas, the correlations are much weaker in the triangular array, 
leading to « 0.3 for high frequencies, and the absence of correlations at low frequencies. 



VI. CONCLUSIONS 



In our numerical study of the conducting nanoparticle films on substrates, we pointed out several steps that play 
an important role in designing the conduction properties of the films. Specifically, we study in quantitative details 
the emergent properties of such nanoparticle assemblies along the line from the self-assembly processes, the topology 
of structures, and conduction through such structures. 

• Self- Assembly processes are leading to different types of structures of the film. We have shown that both regular 
and highly inhomogeneous nanoparticle films, which are typically fabricated in entirely different techniques 0, 
[l9l |. can be assembled within the same process with the bio- recognition binding and varying a single control 
parameter. 

• Conducting paths through the film have been identified and their topological centrality and the dynamical 
flow computed for for V > Vt- Coalescence of flow along different paths contributes to the nonlinear I{V) 
dependences in a range of voltages above the transition Vt- In our simulations, the nonlinearity exponent C 
de pen ds strongly on the structure, in agreement with measured / V curves in different nanoparticle fllms 



• Collective dynamical effects in the tunneling processes deep inside the sample lead to the non-Gaussian current 
fluctuations at the electrode. We have determined various statistical measures of these collective effects. In 
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particular, we find the long-range correlations in charge fluctuations at all nanoparticles, broad distribution 
of flow on junctions through the nano-network, and non-Gaussian distribution of the current increments on 
the electrode. A non-Gaussian distribution have been found in the experiments with nanowires [31]. Our 
numerical results suggest that q — Gaussian distribution with different parameters could be measured in the 
nanoparticle films. The occurrence of the g— Gaussian in a larger class of complex dynamical systems with 
long-range interactions has been discussed recently (3^ . 

Furthermore, we demonstrate how different types of disorder affect the conduction in our comparative study of the 
regular nanoparticle array in the presence of local charge disorder, on one side, and the topologically inhomogeneous 
nano-network without charge disorder, on the other. We hope that our detailed numerical study of the processes 
leading to different structures and the structure-dependent function of the nanoparticle assemblies sheds a new light 
on the problem of engineering of the functional materials with desired characteristics. 
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